A fluctuation relation for weakly ergodic aging systems 
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A fluctuation relation for aging systems is introduced, and verified by extensive numerical simula- 
tions. It is based on the hypothesis of partial equilibration over phase space regions in a scenario 
of entropy-driven relaxation. The relation provides a simple alternative method, amenable of ex- 
perimental implementation, to measure replica symmetry breaking parameters in aging systems. 
The connection with the effective temperatures obtained from the fluctuation-dissipation theorem 
is discussed. 

PACS numbers: 05.40.-a,05.70.Ln 



Non-equilibrium systems are characterized by a net en- 
ergy transfer (in the form of work, heat, or mass) to the 
environment. Aging systems pertain to the category of 
weakly ergodic non-equilibrium systems [l|, Q exhibiting 
slow relaxational dynamics and strong history dependent 
effects where the fluctuation-dissipation theorem (FDT) 
is violated • This has led to the introduction of the 
concept of non-equilibrium or effective temperature 0- 
Despite of the insight gained from exactly solvable models 
and other conceptual attempts (e.g. we still lack a 

clear understanding of the general picture describing ag- 
ing systems. In contrast to stationary systems, aging sys- 
tems are described by two timescales: the waiting time t w 
elapsed since the system was set in the non-equilibrium 
state and the time t > t w at which measurements are 
taken. A characterization of the full spectrum of fluctua- 
tions appears key for a satisfactory understanding of the 
aging state. 

Over the past years several results about energy fluctu- 
ations in non-equilibrium states have been obtained un 
der the heading of fluctuation theorems (FTs) [HI 11 



15(. The new relation is further supported by extensive 
numerical simulations. The existence of an AFR was al- 



FTs take slightly different mathematical forms depend- 
ing on the specific non-equilibrium context [T2, EH . How- 
ever all them share the same common feature: they relate 
probabilities of absorbing and releasing a given amount of 
energy under non-equilibrium conditions; they are useful 
in small systems and short times where energy fluctua- 
tions can be directly measured allowing us, for example, 
to extract free energies of kinetic molecular states from 
irreversible pulling experiments fli| . 

In this work we present a theoretical derivation of a 
fluctuation relation in aging systems. The aging fluctu- 
ation relation (AFR) is based on the hypothesis of par- 
tial equilibration in a scenario of entropy-driven slow re- 
laxation. It can be written in terms of a phase-space 
contraction factor x that bears resemblance to the or- 
der parameter x(q) defined in the context of spin glasses 



ready suggested in Refs. [l6, 17 1 and recently hinted at 
in a quench experiment of a gelatin droplet that exhibits 
a sol-gel transition [l^|. In that reference heat distribu- 
tions were measured and shown to satisfy a fluctuation 
relation for a system in contact with two baths at dif- 
ferent temperatures. Yet it remains unclear whether the 
same relation applies to other aging systems. In con- 
trast, the new relation we are proposing should be gen- 
erally valid in aging systems and is amenable to future 
experimental verification in noise measurements of glass 
formers, critical systems and small systems (e.g. single 
molecules) . 

The aging fluctuation relation (AFR). Consider a sys- 
tem quenched at time t = from an high-temperature 
equilibrium state down to temperature T where it ages 
exhibiting slow relaxation and activated dynamics phe- 
nomena. Such aging state is characterized by low 
entropy-production rate and loss of time-translational in- 
variance [Hj]. During the aging process the system con- 
tinuously exchanges energy with the bath, but some re- 
laxational events result in larger than typical amounts of 
heat Q released to the bath, leading to a net positive en- 
tropy production (AS*) = {Q)/T > 0, where the brackets 
((•••)} denote the average over dynamical histories. 

To analyze the spectrum of heat or entropy production 
fluctuations in the aging state we consider the Bochkov- 
Kuzovlev work fluctuation relation, originally introduced 
as a generalization of the FDT [20(. After a time t w , 
elapsed since the system was quenched, a constant exter- 
nal perturbation of strength h coupled to an observable 
A is applied to the system. This will cause a change 
in A during its subsequent evolution. The entropy pro- 
duction during the time interval [t w ,t] (t > t w ) caused 
by the perturbation is equal to A.St w ,t = Qt w ,t/T = 
h[A(t) - A(t w )]/T = hAA tw , t /T. This quantity has been 
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termed exclusive work in 2l| and satisfies a fluctuation 
relation if the system is equilibrated at t w (which is not 
the case here). AS tmy t is a fluctuating quantity changing 
upon repetition of the same experiment. Moreover in the 
aging regime, where relaxation dynamics is ruled by the 
complex topological structure of the phase space, made 
of many almost degenerate metastable states, ASt m .t dis- 
plays a strong intermittent behavior. This means that 
for fixed t w and AS the requirement AS tm ,t = AS de- 
fines a very broad interval of times t covering many, well 
separated, timcscales. In such a context a data analysis 
for fixed t w and t may be ambiguous since it may mix up 
processes with different timcscales, as noticed some years 
ago in the framework of turbulence 22, HI • To overcome 
this problem we define Pt w {AS) as the probability of ob- 
serving the value ASV^t = AS after t w . 

In Ref. fl6| a similar approach, based on the concept 
of Inherent Structure, was used. The idea was to look 
for ASt mt t associated to the first jump out of an Inherent 
Structure. While this definition could in principle be 
used in a numerical simulation, it is far less useful in a 
real experiment. In contrast, the one proposed here is 
well suited for real experiments. 

Despite the fact that the perturbing field h favors pos- 
itive values of AS tnl , t , trajectories with negative values 
can be observed as well. We argue that, for long enough 
t w , the probability of observing positive and negative 
values of AS satisfies the following fluctuation relation 
(AFR) 



log 
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where is the Boltzmann constant (equal to 1 in the 
following). If at t w the system is in equilibrium, then 
x tw = 1 and Eq. (Q]) reduces to the Bochkov-Kuzovlev 
work fluctuation relation 2(J 21 1 . 

In aging systems, such as glasses, the FDT is not vio- 
lated for times t — t w < t w after switching on the pertur- 
bation. This is because the degrees of freedom whose 
characteristic relaxation times are sufficiently smaller 
than t w have equilibrated. In the present context this 
means that small enough energy transfers between the 
system and the bath occur in equilibrium, and hence 
x tm ~ 1. This remains true as long as \AS\ is smaller 
than a cross-over value AS*, which sets the scale of the 
typical minimum energy transfer in processes that involve 
non-equilibrated degrees of freedom. For \AS\ > AS* 
non-equilibrium energy exchange processes responsible of 
slow relaxation come into play, and the equality xt w = 1 
is violated. Not all degrees of freedom can contribute to 
the relaxation process, some of them being frozen at t w , 
and hence Xt w < 1. As the system ages more degrees of 
freedom equilibrate at t w , and hence AS"* increases with 
t w . In the limit t w 3> t cqi where t oq is the equilibration 
time, all degrees of freedom have equilibrated and Xt w 
converges to 1 for all AS. We note that this may not be 



the case for mean-field models where metastable states 
have infinite lifetime. 

To justify Eq. |T]) we consider a system whose time 
evolution is ruled by the Langcvin equation 



df = —H(ip) 

dt Sip 



(2) 



where ip is an A-dimensional field, £(i) a Gaussian white 
noise (thermal noise) of zero average and correlation 
(£(*)£(*')) = 2TS{t - t'), -6H(<p)/5<p the force aris- 
ing from the conservative energy H(ip), and h the con- 
stant external field coupled to the macroscopic observ- 
able ipt = Si ft ■• where i denotes a site index. Fluc- 
tuation relations derive from the behavior of the prob- 
ability V of a trajectory {f s }s£[t w ,i\ and its reverse 
{fs}s&[t w ,t] = {ft+t w -s}a&[t w ,ty These can be e asily 
computed using the path integral formalism, see e.g. [2JJ: 



n{fs} a& [t w ,t]} = v[{0 s } se[tw , t] \ (3) 



where AS W = ph(ip t - i>t w ) and AS * = f3{H(tp t ) - 
H{ipt w )] = -ln[P c( i(ipt)/P cq {iptJ}- The quantities in 
the exponent depend on the trajectory end-points only, 
then summing over all trajectories from ipt w at t w to ipt 
at time t, and including normalized probability distribu- 
tions Po(pt^) and Pi(ft w ) = Pi (ft) for the initial and 
final states, we get 

e~ AS ™ P(tp t ,t\ip tw ,t w )P (tpt w ) = P(vt w ,t\ft,tw)Pi(ft) 

(4) 

where P((pt,t\ipt m ,t w ) and P(tpt w ,t\cpt,t w ) are the con- 
ditional probabilities of the forward ipt m — > tpt and re- 
verse ipt —> ipt w trajectories, respectively, and AStot = 
AS tvj , t -AS c * + AS hl with AS b = -ki[P 1 (<p t )/P (<p t J]. 
From this relation the following identity follows 



P W (-AS) = (S(AS tw , t ~ AS) e~ As ^) 



t m ,t 



(5) 



where the average is over the forward process tp tw — > ipt 
with initial probability distribution Po(tp tnl ). 

After the quench the system partially equilibrates in- 
side independent phase-space regions (that we will call 
cages), from which it will escape only after a time t—t w ~ 
t w . Therefore, when t—t w <^t w the system is in (partial) 
equilibrium with the thermal bath, Po(<p) = Pi(f) oc 
P cq (<p) so that AS cq = AS b , AStot = AS tw ,t, and from 
© one gets eq. (fTJ) with x tul = 1. 

To study the opposite limit t — t w 3> t w , where the 
system can access different cages, we observe that ASt w ,t 
depends only on the macroscopic variables tpt and ipt„ , 
then the average on the r.h.s of ([5|) can be done by partial 
classification, that is by averaging first over all paths with 
given initial and final states, and then over ip t and tpt w '■ 



P twtt (-AS) = J dip tw J dip t 5(AS twtt -AS) 



x ( e -AStot) 



(6) 
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where ((• ■ -))ii> tw t m ;i)tt denote dynamical averages re- 
stricted to those trajectories starting with ipt m at t w 
and ending with ip t at t. Assuming the system is par- 
tially equilibrated over cages, the probability P(ip\tp) 
of a state tp in a cage with fixed ip is proportional to 
P eq (ip) times the probability of having -0 in a cage. Thus 
P(y#) oc P oq (^) x n cagc (il))/n tuj {ip), where QcageWO is 
the number of states with ip inside the cage, divided by 
the total number Qt w (tp) of accessible states, not neces- 
sarily in the same cage, with ip. Under the hypothesis 
of partial equilibrium S , cage ('0) = hiJ7 cage (0) is the ther- 
mal equilibrium entropy in the cage. Q, tm (tp) depends 
on system age since more we wait more degrees of free- 
dom relax. S tm (tp) = In Q, tw (tp) is then smaller than the 
full thermodynamic entropy and converges to it only for 
t w 3> icq- Using this ansatz for the PDF of the initial 
and final states, and the relation dS cage (tp)/dtp = (3h 
together with the analogous dS tm (tp)/dtp = xt w f)h, cor- 
rected through the coefficient xt w < 1 to account for 
the frozen degrees of freedom, we have AS cq — A 5b = 
/3h(tp t -tpt a )-xt w /3h(tp t -tpt w ) = (l-x t JASt w ,t- Insert- 
ing this form into (JB|) the AFR JTJ) follows. The coefficient 
Xt w measures the phase space contraction due to frozen 
degrees of freedom at t w . Clearly x tvj —> 1 as t w » t cq . 

Numerical tests. We have tested the AFR ([T]) through 
Monte Carlo simulations of several model systems, but 
we report results for only three of them. The system, ini- 
tially prepared in an high temperature equilibrium state, 
is instantaneously quenched to a temperature T below 
the freezing transition temperature. After t w a pertur- 
bation of small intensity, to ensure a good statistics for 
trajectories with AS < 0, is applied and the fluctua- 
tions AA tmt t of the conjugated variable A are recorded 
at fixed time intervals t — t w . The maximum recording 
time t was taken much larger than t w to ensure good 
statistics. The procedure was repeated several times, 
P tm (AS) calculated from data binning and Eq.([T]) tested 
to extract the value of Xt w ■ To compare it with the 
parameter x = T/T c ff(t w ) obtained from the FDT, the 
fluctuation-dissipation (FD) plots in time-domain were 
also measured using the standard procedures Q ■ Both x 
parameters (the one derived from the AFR, Eq.([T]), and 
the one derived from FDT) asymptotically coincide un- 
der general assumptions, see [25[. 

The first model is the Random Orthogonal model 
(ROM) defined by the Hamiltonian [1] 



H 



E 

l<i<j<N 



J ij &i &j i 



(7) 



where <Tj = ±1 are Ising spins and Jij = Jji quenched 
Gaussian variables of zero mean and variance 1/N sat- 
isfying J2k^ikJkj = lQSij, with Ju = 0. The system 
is perturbed by a uniform magnetic field of strength h 
conjugated to the total magnetization M(t) = J2i a i(t)- 
This model describes structural glasses in the mode 
coupling theory (MCT) approximation and has a dy- 



namical MCT transition at Td = 0.536. Below Td the 
system is dynamically confined into one of the many (ex- 
ponentially large in number) metastable states and can- 
not reach full equilibrium. The equilibrium transition 
occurs at the lower static (or Kauzmann) temperature 
T c = 0.25. The low-temperature behavior of the model is 
described by a one-step replica symmetry (1RSB) break- 
ing order parameter. Typical signature of this is a two- 
slopes FD plot. 
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FIG. 1: Numerical test of the AFR in the ROM. Main 
plot is a test of Eq.Q for the model with N = 1000, T = 0.2 
and h = 0.1 at three values of t w . The dashed line corresponds 
to x tw = 1 while full lines to x tw = 0.271,0.287,0.299 for 
t w = 64, 1024, 16384 respectively, (a) P iw (AS) and (b) FD 
plot. The dashed and continuous lines shown in the FD plot 
have slopes equal to those shown in the main plot, (c) Zoom 
of the region corresponding to xt w = 1 (intra-cage relaxation). 

Figure [T] shows results for the ROM. Two regimes can 
be distinguished, x tw = 1 for |A5| < AS* and x tw < 1 
for \AS\ > AS* (AS* ~ 2 for t w = 1024). The values of 
Xt m agree quite well with x — T/T e g, extracted from the 
FD plot [inset (b)]. 

As a more realistic system we have studied a 80 : 20 
binary mixture of type A and B particles interacting via 
a Lennard- Jones pair potential (BMLJ): 



V a p(r) = 4e Q/3 



(tT-G?) 



(8) 



where a, /3 = A,B, r is the distance between the two 
particles and the parameters a a p , e a p stand for the ef- 
fective diameters and well depths between species a, (3. 
The parameters for length and energy measured in units 
of oaa and eAA arc £bb = 0.5, 6ab = 1.5, <tbb = 0.88 
and cab = 0.80, and are taken to prevent crystallization 
[III. With this choice a system of reduced density p = 1.2 
exhibits a glass transition well described by the MCT at 
the critical temperature Tmct — 0.435. The study of the 
AFR was done by adding at time t w an external poten- 
tial of the form U ■ €j cos(fc • r), where Vo < T and ej 
arc i.i.d. (quenched) random variables equal to ±1 with 
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equal probability, and recording the conjugated observ- 
able Ak{t) = ^2j6jexp[ik ■ rj(t)] ■ Results for the 
AFR arc shown, and compared with standard FD plots, 
in Fig [2] Also in this case the agreement between the 
x extracted from AFR and that from FD plot is rather 
good. Interestingly curves in the AFR plot do not ex- 
hibit the bending that can be seen in FD plots. Such 
bending is due to finite-size effects that for short t w pro- 
duce an equilibration of the system, and a depletion of 
the statistics of rare events for longer t w . The plot then 
bends upward for short t w and downward for longer t w . 




FIG. 2: Numerical test of the AFR in the BMLJ 
model. Main plot is a test of Eq.(ff]) for the model with 
parameters N = 500, Vb = 0.1, |fe| = 7.25, T = 0.3 at 
three values of t w . The dashed line corresponds to xt m = 1 
while full lines (red, green, blue) to x± w = 0.40, 0.54, 0.60 for 
t w = 10 2 , 10 3 , 10 4 respectively, (a) Pt w (AS) and (b) FD plot. 
The dashed and continuous lines (red, green, blue) in the FD 
plot have slopes equal to those shown in the main plot. 

Both above systems are described by a two-steps relax- 
ation, or 1RSB, scenario. As last example we have con- 
sidered the 3-dimensional ± J Edwards- Anderson model 
(±J-EA) defined by the Hamiltonian 0, but with Jy 
randomly chosen equal to ±1 if the sites j and i are 
nearest-neighbors on a cubic 3-dimensional lattice, and 
zero otherwise. Numerical investigation indicates that 
below T ~ 1.14 there is a spin-glass phase described 
by a continuous-step relaxation, or Full-RSB, scenario. 
Typical signature of this is a continuous-slope FD plot. 
Numerical tests of Eq.([T]) for the ±J-EA at two differ- 
ent t w are shown in Fig 02 together with the FD plot. 
The AFR is well verified also in this model. However, as 
shown in inset (b) and in contrast to the previous 1RSB 
models, the phase contraction coefficient Xt w agrees with 
T/T e s only at the early stage, where the FD plot de- 
part from equilibrium. This is in agreement with the 
aforementioned argument, according to which Xt w gives 
information on the phase-space partition at t w . 

Discussion. Summarizing, the AFR (UJ shows a 
promising route to experimentally test the partial 




FIG. 3: Numerical test of the AFR in the 3D ±J EA 

model. Main plot is a test of Eq.fTJ for the model with 
parameters L = 16, T = 0.7 and h = 0.1 at two values of 
t w . The dashed line corresponds to x± w = 1 while full lines 
(red,blue) to x tw = 0.383, 0.450 for t w = 95, 1020 respectively. 
(a) Pt w (AS) and (b) FD plot. The dashed and continuous 
lines (red, blue) in the FD plot have slopes equal to those 
shown in the main plot. 



equilibrium- entropy driven scenario in slowly relaxing 
systems from noise measurements. A theoretical deriva- 
tion of the AFR was given and its validity verified by ex- 
tended numerical experiments. The connection between 
the values of xt w extracted from the AFR and FD plots 
was shown. Remarkably enough, and in contrast to FD 
plots, extracting the value of Xt w does not require mea- 
suring aging correlation functions. By only measuring 
the statistics of AS" to an externally applied perturba- 
tion, P tw (AS), we can test the validity of Eq.([T]) to ex- 
tract the value of Xt w ■ We emphasize that in order to 
test Eq.([T]) statistical events with AS < must be ob- 
served. Since the average value of AS continuously in- 
creases with t, only rare events with AS < give full 
meaning to the AFR. A similar situation is encountered 
in the Gallavotti- Cohen theorem for steady state sys- 
tems [3(J. Eq.([T|) is ready to be employed in mesoscopic 
systems (e.g. magneto-conductance fluctuations in spin 
glasses and electron glasses [3l|, |32| ) and single molecule 
experiments. The latter include molecular systems ex- 
hibiting slow folding due to disorder and frustration (e.g. 
RNA) or slow binding kinetics (e.g. peptides or proteins 
binding DNA). Ultimately, small systems may provide a 
direct access to experimentally measure the always elu- 
sive spin-glass order parameter. 
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